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ABSTRACT 


We study the influence of the presence of a strong bar in disc galaxies which 
host an active galactic nucleus (AGN). Using data from the Sloan Digital Sky Survey 
and morphological classihcations from the Galaxy Zoo 2 project, we create a volume- 
limited sample of 19,756 disc galaxies at 0.01 < z < 0.05 which have been visually 
examined for the presence of a bar. Within this sample, AGN host galaxies have a 
higher overall percentage of bars (51.8%) than inactive galaxies exhibiting central star 
formation (37.1%). This difference is primarily due to known effects; that the presence 
of both AGN and galactic bars is strongly correlated with both the stellar mass and 
integrated colour of the host galaxy. We control for this effect by examining the dif¬ 
ference in AGN fraction between barred and unbarred galaxies in hxed bins of mass 
and colour. Once this effect is accounted for, there remains a small but statistically 
signihcant increase that represents 16% of the average barred AGN fraction. Using the 
A[o iii]/AfBH ratio as a measure of AGN strength, we show that barred AGN do not 
exhibit stronger accretion than unbarred AGN at a fixed mass and colour. The data are 
consistent with a model in which bar-driven fueling does contribute to the probability 
of an actively growing black hole, but in which other dynamical mechanisms must 
contribute to the direct AGN fueling via smaller, non-axisymmetric perturbations. 


1 INTRODUCTION 


Supermassive black holes exist at the centres of most (if not 


all) massive galaxies ( 

Kormendy & Richstone 

1995 Rich- 

stone et al.||19981 |Kormendy & Gebhardt||2001 

|Ghez et al. 

2008 

). The evolution of the black hole is closely tied to that 


of the host galaxy; hence, understanding the conditions that 
drive black hole growth is key for a complete picture of 
galactic evolution. While most black holes are not actively 
growing, a small fraction are observed to accrete matter 
and cause the surrounding material to emit powerful pan¬ 
chromatic radiation. The central region of a galaxy which 
encompasses these “active” black holes, along with the sur¬ 
rounding accretion disk and ionized gas clouds, is an active 
galactic nucleus (AGN). Since the bolometric luminosity of 
the AGN can be comparable to (or greater than) the in¬ 
tegrated stellar luminosity (as high as L ~ 10"^^ erg s“^) 


the black holes have an important effect on the host galaxy, 
controlling the amount of star formation via AGN feedback, 
as well as contributing toward the net reionization of the 
intergalactic medium ( Heckman fc Best|2014 |. Understand¬ 
ing the fueling mechanism(s) for AGN is thus critical for 
studying galaxies, both in the nearby Universe and at higher 
redshifts. 


The precise physics that govern the relationship be¬ 
tween AGN and their host galaxies is an area of intense 
study. This includes the AGN fueling mechanism — while 
there is strong evidence that there is sufficient gas in the ISM 
to keep the accretion disc supplied with enough material to 
radiate at typical bolometric AGN luminosities ( [Shlosm^ 
et al. JT^I 1990[ |, the dynamical mechanisms that drive the 
gas within the black hole’s sphere of influence are difficult 
to observe directly, especially at extragalactic distances. In 
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order to initiate (or continue) AGN activity, gas must lose 
enough angular momentum in a short timeframe to reduce 
its orbit from scales of kiloparsecs down to parsecs. |Shlos-| 


man et al. (19891 analytically showed that while gas can 


lose angular momentum due to turbulent viscous processes, 
these are too slow to be the only mechanism involved. Later 
N-body simulations have shown viscous torques on the gas 


are negligible and do not directly initiate inflows (Bournaud 


et al. 20051, further arguing for an additional method of 


radial gas transport. 

One possibility is that the presence of a large-scale bar 
may supplement viscous torques and further drive AGN fuel¬ 
ing. Bars efficiently transport angular momentum within the 


disc (Athanassoula 2003 Kormendy & Kennicutt 20041, and 


are ubiquitous features in disc galaxies in the local Universe 


(|Eskridge et al. 

2000 [Laurikainen et al. 2004 

Menendez- 

Delmestre et al. 

|2007| Masters et al. 

[2011 G 

leung et al. 


2013). Simulations ( Athanassoula|19^ Friedli & Benz 1993 

Ann fc Thakur||2005[ ) show that stellar bars, whose lengths 

are on the order of kiloparsecs, do drive gas into the circum- 
nuclear region (scales of 100 pc) of galaxies; observational 
studies have also shown an increase in the amount of central 
star formation for barred galaxies ( Ellison et al.|20TT ). This 
combination of simple analytical models, simulations, and 
observations clearly points toward galactic bars preferen¬ 
tially driving gas to the centres of their galaxies. It is still an 
open question, though, whether this gas is ultimately driven 
to the central 1 — 10 pc scales, which theoretical models sug¬ 
gest are required for accretion around the central black hole 
of the AGN. 

Theoretical models for alternate modes to bar-driven 
fueling also exist. Numerical simulations from [Hopkins fc| 
IQuataert] ( |2010[ ) examine several possible mechanisms be¬ 
hind angular momentum transport for a range of galaxy 
morphologies (bars, spirals, rings, clumpy and irregular 
shapes, mergers) at different radial scales. For each morpho¬ 
logical type, gas transported from larger to smaller (~ 1 kpc) 
radii “piles up” due to decreasing efficiency in the processes 
that induce torque. If this pile-up of gas is sufficiently mas¬ 
sive, it becomes self-gravitating and can efficiently transport 
angular momentum down to scales of ~ 10 pc. This “stuff 
within stuff” model is similar to the second half of IShlosmanI 


et al. (19891’s “bars within bars” model. The difference is 


that the “bars within bars” model assumes that a large- 
scale bar is the primary mechanism that transports the gas 
inward to form the gaseous disc, while [Hopkins fc Quataert] 
(20101 show that many large-scale morphologies are capable 


of producing a secondary instability and fueling an AGN, 
suggesting that this process may not be restricted to classic 
large-scale bars. 

Many studies have focused on observational correlations 
between the presence of a galactic bar (typically identified 
at optical wavelengths) and that of an AGN (identified by 
optical line ratios or widths). Some studies (eg. Ho et al. 


1997 Mulchaey &i Regan||1997| Hunt fc Malkan||i999| ) find 

similar bar fractions for both AGN and inactive galaxies 
and hence report no correlation. The significance of these 
fractions, however, is hindered by small sample sizes, typi¬ 
cally with fewer than 100 barred AGN hosts. More recent 
studies (Knapen et al.pOOO Laine et al. 2002 Laurikainen 


et al. 2004) report increases of 20 — 23% in the bar frac¬ 


tions for AGN when compared to non-AGN hosts. Despite 


larger numbers of AGN, the results are still only significant 
at the 2.5(7 level. Rather than comparing the likelihood of 
active and inactive galaxies to host bars, as is most common 
among previous studies, Cisternas et al. (2013) accounted 
for a continuum of values by quantifying bar strength and 
activity level in local X-ray identified AGN. While no corre¬ 
lation was found, these data probe only the low-luminosity 
AGN regime (Lx ~ 4 x 10^* erg s“^). In the high redshift 
universe, Gheung et al. (2015) find no compelling evidence 


that bars are more likely to lie in AGN hosts than non-AGN 
hosts. 


Several recent studies have focused on optical identi¬ 
fications of bars and AGN, primarily using data from the 
Sloan Digital Sky Survey (SDSS). We compare these meth¬ 
ods and results in Table [T| Among these studies, neither [Leej 


et al. (2012) nor Martini et al. (2003) find any correlation 


between the presence of strong galactic bars and AGN, but 
do not rule out the possibility of smaller, nuclear bars in¬ 


fluencing AGN activity. In contrast. Oh et al. (2012); Hao 


et al. (20091; Alonso et al. (20131 all find evidence of bar ef¬ 


fects in AGN — however, they disagree on both the strength 
of the effect and whether it affects both black hole fueling 
and/or central star formation. One possible reason for the 
discrepancy is the lack of a consistent scheme for classifying 
AGN. While the BPT diagram based on optical line ratios 
( Baldwin et al.|fl981| is among the most common methods 
for identifying AGN, the demarcation between star-forming 
and AGN host galaxies is not consistent; some use thejKew^ 


ley et al. (2001) criterion that excludes composite galaxies, 


while others use Kauffmann et al. (2003a) and include these 
along with Seyferts as AGN. The inclusion of LINERs can 
also complicate the picture; the high line ratios in at least 
some LINERs are spatially extended and thus likely of a 


non-AGN origin (Sarzi et al. 2010 Yan & Blanton 2012 
Singh et al.|2013 l. 


Other challenges result from the task of identifying 
galactic bars, which is often done by visual inspection of 
optical images by individuals or small groups of experts. 
This introduces potential complications when there is dis¬ 
agreement between classifiers, especially in the cases of weak 
or nuclear bars. With only a single (or a few) classifications 
per image, such disagreements are difficult to resolve. Fur¬ 
thermore, individual visual inspection can limit the effective 
sample size due to the amount of time required to inspect 
images one by one. Our work avoids these problems by using 
crowdsourced citizen science classifications to identify galac¬ 
tic bars, where many individuals (an average of 27 classifiers 
for bar detection in this study) analyze each galaxy, and the 
presence of a bar is quantified as a calibrated vote fraction. 


This paper re-examines the relationship between bars 
and AGN in disc galaxies by using Galaxy Zoo morpholog¬ 
ical classifications, and by using a strict AGN classification 
scheme which only selects Seyfert galaxies. We use this data 
to consider three physical scenarios for describing the role 
bars may (or may not) play in AGN fueling: I) Bars are 
necessary to fuel AGN, H) Bars are one of several ways to 
fuel AGN, or HI) Bars do not fuel AGN. We discuss each 
of these possibilities in Section and suggest the means by 
which the existence barred AGN, unbarred AGN, barred 
non-AGN, and unbarred non-AGN may be explained within 
the context of each model. We then report the scenario which 
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• 

pbar=0.00 

58774170834360243 

Pbar=0.11 

58772986068639824 

Pbar=0.25 

58772892417241196 






• 

Pbar=0.35 

587726016161448090 

Pbar=0.43 

588297863115374695 

PbBr=0.55 

587738615414718507 

s 



• 



pbar— 0.65 
587736586050207931 

pbar=0.79 

587724648184938507 

pbar=1 .00 

587728669874126911 


Figure 1. Examples of the SDSS images used in Galaxy Zoo 2, 
sorted by increasing Pbar (the weighted percentage of users that 
detected a bar in each image). All galaxies are from our final 
analysis sample of “not edge-on” disc galaxies. The white lines in 
the upper left of each image represent a physical scale of 5 kpc. We 
also give Pbar ^^nd the SDSS objectIDs for each galaxy. Top row: 
Galaxies with Pbar< 0-3, which in this paper are designated as 
unbarred. Middle and bottom rows: Galaxies withpbari? 0-3, 
which we designate as reliably barred. 


we find to be best supported by both our observations and 
current theoretical models and simulations. 

In Sectionj^we describe our sample selection. Sectionj^ 
includes our data, with mass and colour distributions of the 
different activity types, both barred and unbarred, as well 
as a comparison between accretion strengths of barred and 
unbarred AGN. Interpretations of these results are discussed 
in Section]^ and the main conclusions are outlined in Sec¬ 
tion We adopt a ACDM cosmology throu ghout the paper 
of Qm = 0.27 and Ho = 71 km s“^ Mpc“^ (Planck Collab- 
oration et al.|[2M3|. 


2 DATA AND SAMPLE SELECTION 


Our parent sample of galaxies is taken from the SDSS 
Data Release 7 ( Abazajian et alT]|2009 l. From the spectro¬ 
scopic Main Galaxy Sample (Strauss et al. 20021, we se¬ 
lect galaxies within the redshift interval 0.01 < z < 0.05 
— the lower limit excludes galaxies whose angular size sig¬ 
nificantly exceeds the spectroscopic fiber, and the upper 
limit is chosen so that a reasonable estimate of bar detec¬ 
tion can be made by visual inspection. From this, we create 
a volume-limited sample by applying an additional cut of 
,petro —19.5 AB mag. 

Within the volume-limited sample, we use morpholog¬ 
ical cuts to select only disc galaxies at low inclination an¬ 
gles that are candidates for the presence of galactic bars 


(described below). These cuts result in the final sample of 
19,756 disc galaxies used in the remainder of this paper. 


2.1 Bar classifications and Galaxy Zoo 2 


To select disc galaxies and measure the presence of a 
bar, we use data from the online citizen science project 
Galaxy Zoo 2 (GZ2)[^With the help of over 80,000 volun¬ 
teers providing over 16 million classifications of over 300,000 
galaxies, Galaxy Zoo 2 is the largest extant survey of de¬ 
tailed galaxy morphology. Volunteers are shown colour im¬ 
ages of galaxies taken from the SDSS (Figure [^, and are 
then prompted through a decision tree in which they an¬ 
swer questions about the galaxy’s structure. For a detailed 
discussion on the Galaxy Zoo 2 project and its decision tree, 
see 


Willett et al. (20131. 


Since bars only appear in disc galaxies, the sample must 
be limited to disc galaxies in which a bar can be seen via 
visual inspection. We begin by selecting galaxies for which 
at least 10 people answered the question, “Is there a sign of 
a bar feature through the centre of the galaxy?”, thus reject¬ 
ing vote fractions with low statistical significance. Because 
questions in GZ2 are implemented as part of a decision tree 
( [Willett et al.|20T^ , users must have identified a galaxy as a 
disc and as not edge-on before answering the bar question. In 
this way, the cut of Nbar (S 10 increases the likelihood that 
the galaxy in question is a candidate for having a bar. This 
cut is not complete, however, for galaxies which have a high 
number of total classifications. In these cases, the number of 
users to answer the bar fraction may still be small compared 
to the number of users identifying the galaxy as either not 
disc-like, or as an edge-on galaxy. Therefore cuts are also 
applied to the vote fractions relating to questions preceding 
the bar question. The first question of the GZ2 tree reads, 
“Is the galaxy simply smooth and rounded, with no sign of 
a disc?” [Willett et al.| ( [2013[ ) determined the threshold frac¬ 
tion of “features or disc” answers required to classify the 
galaxy as a disc, when combined with the cut Nbar -S 10, 
to be ^features or disk^ 0.227. We emphasize that the cuts 


provided in Willett et al. (20131 are intended to be mini¬ 


mum values for determining well sampled galaxies. We thus 
chose to adopt a slightly higher value ofpfeatures or disk^ 0.35 
to create the cleanest possible sample, based on a visual in¬ 
spection of a subsample of galaxies with these cuts. To assess 
whether the results would be affected by this choice, we also 


created a sample with the original Willett et al. (20131 cuts. 


This choice increased the number of AGN in the sample by 
24, and did not affect the final results. Therefore we present 
the sample using our more conservative cuts in this paper. 

Following an answer of “features or disc” for the first 
question, the volunteer is then asked “Could this be a disc 
viewed edge-on?” Bars become increasingly difficult to de¬ 
tect in galaxies at high inclination angles, and are nearly im¬ 
possible to detect in edge-on galaxies without careful isopho- 
tal mapping. The threshold vote fraction determined by 
Willett et al. (20131 of a “No” answer to this question is 


Pnot edge-on^ 0.519. We again adopt a slightly more conser¬ 
vative value of Pnot edge-on^ 0.6 based on visual inspection 
of a subsample. The combination of feature/disc galaxies 


^ zoo2.galaxyzoo.org 
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Martini et al. ( 2003 1 

Hao et al. ( 2009 1 

Lee et al. 12012 l 

Oh et al. (2012 1 

Alonso et al. (2013 1 

This work 

Redshift range 

2 : < 0.038 

0.01 < 2 < 0.03 

0.02 < 2 < 0.055 

0.01 < 2 < 0.05 

2 < 0.1 

0.01 < 2 < 0.05 

Abs. magnitude 

range 

Bt < 13.4 

18.5 < Mg < 
-22.0 

Mr < -19.5 ± 
51og(6) 

Mr < -19 

Mg < -16.5 

^z,petro 19.5 

Inclination limit 

R 25 < 0.35 

i < 60° 

h/a > 0.6 

b/a > 0.7 

h/a > 0.4 

Pnot edge —on^ 6.6 

AGN classification 
method 

varied 

FWHM(Ho) > 
1200 km/s and 
Ka03 

KeOl 

Ka03 

Ka03 

S07, WISE 

AGN type(s) 

Type 1 and 2 
Seyferts, LINERs 

Type 2 Seyfert, 
LINER, compos¬ 
ite 

Type 2 Seyfert, 
LINER 

Type 2 Seyfert, 
LINER, composite 

Type 2 Seyfert, 
LINER, composite 

Type 2 Seyfert 

Bar classification 

method 

visual inspection 

ellipse fitting 

visual inspection 

visual inspection 

visual inspection 

crowdsourced 
visual inspection 

Number of AGN in 
sample 

28 

128 

1742 

1397 

6772 

681 

Fraction of AGN 
hosts that are 

barred 

28.6% 

47% 

49% 

51% 

28.5% 

51.8% 


Table 1. Summary of recent studies comparing the presence of galactic bars and active galactic nuclei, including new results from this work. 
[Martini et al.| < |2003| l is the only study with neither uniform selection criteria for galaxies nor a volume-limited sample. AGN classifications from 
optical line ratios and the BPT diagram are separated by the following demarcations: KeOl = [Kewley et al!] l |2001[ |; Ka03 = [Kauffmann et al.j 
< |2003a| l; S07 = |Schawinski et al.| l |2007| . 


that are not edge-on for these two thresholds results in the 
final sample size of 19,756 galaxies used in this paper (Ta- 
ble|^. 

As a check that our selection of “not edge-on” disc 
galaxies can be reliably used to identify a bar, we examine 
the inclination angle of the sample, which is approximated 
by the ratio of the best fit of the semi-major and -minor 
axes i = cos” ^ (a/6) as measured in r-band by the SDSS 
pipeline. Figurej^shows the strong correlation between i and 
Pnot edge-on, with a sharp cutoff near i = 70°. Our cutoff of 
Pnot edge-on^ 0.6 effectively limits the sample to inclination 
angles of i < 67°. In Figure we also show the dependence 
of the GZ2 bar fraction on pnot edge-on. The bar fraction re¬ 
mains roughly constant (±0.1) between 0.3 <Pnot edge-on< 
1.0 and drops to zero at pnot edge-on< 0.1. Since the true 
bar fraction is expected to be independent of i (a purely 
geometrical effect assumed to have a random distribution), 
any change in the bar fraction would reflect the ability of vi¬ 
sual inspection to detect a bar in a highly inclined disc. The 
constant bar fraction out to our limit of Pnot edge-onJ? 0.6 
(and well beyond) is a necessary requirement for an unbi¬ 
ased selection of barred galaxies; as a result, we are confident 
that the crowdsourced bar classifications in this sample are 
reliable. 

Finally, if the volunteer answers “No” to the edge- 
on question, they are asked “Is there a sign of a bar fea¬ 
ture through the centre of the galaxy?” Possible answers to 
this question are either “Bar” or “No bar”. [Willett et al.| 
(|2013|) compared expert classifications of barred galaxies 


from both Nair & Abraham (20101 and Baillard et al. (20111 


to Galaxy Zoo 2 data, and show that a threshold of Pbar^ 0.3 
is the most reliable separator of the barred from unbarred 
population (see their Figure 10). We adopt the same thresh¬ 
old of Pbar^ 0.3 for determining whether a galaxy has a bar 
(see Figure for images of galaxies with different values of 

Pbar). 


We compare our morphology cuts to those used byjMas-j 


ters et al. (20111, who used an early release of GZ2 data to 


identify barred galaxies. Their study also required Nbar 1? 10 
and claim that this cut alone is sufficient to restrict the 
sample to disc galaxies without applying an additional cut 
on pfeatures or disk. This assumption was reasonable at the 
time since the Galaxy Zoo 2 project was still collecting 
data, and the number of classifications per galaxy was lower 
than in the final catalog. The median number of classifica¬ 
tions per galaxy is roughly 30% higher, and so our data is 
more susceptible to contamination by non-disc galaxies with 
high classification counts. This makes an additional cut on 
Pfeatures or disknecessary. To remove edge-on discs, jMastersj 


et al. (20111 set an inclination limit of log(a/6) < 0.3, or 
i ~ 60°; this is comparable to our pnot edge-on cut, which 
corresponds to roughly i ~ 67°. To select barred galaxies, 
a majority vote fraction of pbar> 0.5 was required, higher 
than our value ofpbar^ 0.3. We are nevertheless confident in 
our threshold which was determined by the more recent and 
detailed analysis of the GZ2 data by Willett et al. (20131 
as described above. Additionally, the data released at the 
time of [Masters et al.] (|2011| had not yet been reduced via 
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Figure 2. Left: Fraction of “not-edge-on” votes vs. inclination angle {i = cos“^[a/6]) for the disc galaxies in our GZ2 sample. An angle of 
0° means the galaxy is completely face-on, while 90° is completely edge-on. GZ2 users consider a galaxy as “not edge-on” if the inclination 
angle is less than i ~ 70°. Right: Fraction of barred galaxies vs. fraction of “not edge-on” galaxies. The bar fraction is independent of 
the edge-on degree of the galaxies (above pnotedgeon ~ 0.3); the ability of users to detect bars does not decrease with inclination until 
Pnotedgeon 0.3, or i ~ 70°. Error bars are 95% Bayesian binomial confidence intervals ( |CamerQn[2011[ ). This demonstrates that GZ2 
data can reliably identify bars even in moderately-inclined disc galaxies. 


weighting and debiasing; these differences in vote fractions 
also contribute to the different cuts used in our study. 


2.2 Activity type classification 


We use flux measurements from the 2012 release of the 
OSSY catalogue pr et al. MB to classify disc galaxies 
as either star-forming, composite, AGN, LINER, or quies¬ 
cent (also known as “undetermined”). This method employs 
ratios of [O iii]/H/3 fluxes as a function of [N ii], [S ii], or 
[O l] over Ha according to the BPT diagnostics. Our method 
for selecting AGN is the same as used by |Schawinski et al.| 
| |2007[[2C)10| . First, we use the [N ii]/Ha ratio (Figure |^). 
Any galaxy that does not have S/N > 3 for any of the four 
lines is unclassifiable via this method (possibly due to be¬ 
ing gas-poor) and labeled “undetermined.” Next, any galaxy 
which falls below the Kewley et al. (20011 extreme starburst 


line is classified as star-forming, and those that fall between 
this and the [Kauffmann et al)] ( |2003a[ ) empirical starburst 
line are classified as composite. We note that some of these 
composite galaxies may be potential AGN, but we cannot 
cleanly separate the AGN contribution from star formation 


and thus exclude them from our sample (Schawinski et al. 


2010 ). 


Next, we identify the remaining galaxies (above the 
extreme starburst line) as either Seyfert AGNs or LIN- 
ERs. [Kewley et "n? (2006) showed that both [O l]/Ha and 
[S ll]/Ha diagrams are better-suited to distinguish AGN 
from LINERs; we thus use diagram (c) in Figure]^ if these 
galaxies also have S/N > 3 in [O l]. For galaxies which do 
not have S/N > 3 in [O l], but do in [S ll], we use diagram 


(b). In both cases, we use the AGN-LINER division line of 


Kewley et al. (20061. For the remaining galaxies, we use di¬ 


agram (a) and implement the AGN-LINER division line of 


Schawinski et al. (2007). 


Finally, to detect any AGN that may have been op¬ 
tically mis-classified due to obscuration, we identify AGN 
based on their infrared continuum shape using data from 
the Wide-field Infrared Survey Explorer ( Wright et al.|2010{ 
WISE). We identify as an AGN any galaxy with {W1 — 
W2) ^ 0.8 ( |Stern et al.||201^ . Based on infrared data, we 
re-classified fourteen galaxies (originally classified optically 
as three star-forming, ten composites, and one LINER) as 
AGN. 


We show the results of the activity type and morpho¬ 
logical classifications in Table The numbers and fractions 
of each activity type with respect to the full sample are 
shown, as well as the numbers and fractions of barred galax¬ 
ies within each activity type. These results are discussed in 
Section |3 


RESULTS 


To determine whether a correlation exists between galaxies 
that host an AGN and those that contain large-scale stellar 
bars, we examine the fractions of barred and unbarred AGN 
with respect to mass, colour, and AGN strength. We use stel¬ 
lar masses from the AVERAGE values in the MPA-JHU DR7 
catalogue ( Kauffmann et aL]|2003b ). Colours are ° °(w — r) 
values from SDSS DR7, which have been both de-reddened 
for Galactic extinction and fc-corrected to redshift 2 = 0.0 
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-1.0 -0.5 0.0 0.5 -0.8 -0.4 0.0 0.4 -2.0 -1.5 -1.0 -0.5 

log ([Nil] A 6583/Ha) log ([SII] A A 6716,6730/Ha) log ([01] A 6300/Ha) 


Figure 3. Optical line diagnostics for activity types of 19,756 disc galaxies. Any galaxy with S/N < 3 for [O ill], H/3, [N ii], or Ha 
is unclassifiable using this method and labeled as “undetermined”. The 3,619 undetermined galaxies do not appear on the diagram 
above. The remaining 16,137 galaxies were categorized according to the above diagrams in the following order, based on the method of 
[Schawinski et ah] l 2007[ l. First, diagram (a) was used to identify star-forming and composite galaxies. Any galaxy below the Ka03 line was 
classified as star-forming, while those that fell between the Ka03 and KeOl lines were classified as composite. Next, to distinguish AGN 
from LINERs, we use diagrams (b) and (c). If a galaxy had S/N > 3 for [O i], diagram (c) was used. If a galaxy did not have S/N > 3 
for [O i], but did for [S ii], diagram (b) was used. Last, if a galaxy did not have S/N > 3 for [O i] or [S ii], but did for [N ii], diagram 
(a) was used. In each panel, only galaxies with S/N > 3 for all four lines required by that diagram are shown. Galaxies designated AGN 
by any of the three optical line diagnostics are plotted as blue points, while the black shading represents the full sample of emission-line 
galaxies. 


All discs 


Barred discs 


Activity type Number /total(%) 


Number /bar(%) 


barred subsamples. The median AGN is more massive (by 
0.6 dex) and redder (by 0.5 mag) than the median star- 


star-forming 

composite 

AGN 

LINER 

undetermined 


11282 

2853 

681 

1321 

3619 


try -1 -F0.7 
-0.7 

14 4 

-0.4 
o / -hO.3 
- 0.2 
6 7 

-0.4 
+0.6 
-0.5 


4183 

37.1 

AGN and star-forming galaxies in the local Universe ( 

Schaw- 

1301 

45.6 +]® 

inski et al.||2007| |Lee et al.||2012| |Oh et al.||2012| 

Alonso 

353 

51.8 + ® 

et al.|2013l 

Aird et al. ( 

2012 

1 demonstrate that this differ- 


18.3 


695 

1654 


52.6 

45.7 


+2.7 

-2.7 

+ 1.6 

- 1.6 


total 


19756 


100 


8186 


41.4 


+0.7 

-0.7 


Table 2. Results of activity classification for our sample of 19,756 
not edge-on disc galaxies, /total is the percentage of the total 
sample represented by each activity (number of galaxies of that 
type / total number of galaxies), /bar is the percentage of each 
subsample that are barred (number of galaxies of that type that 
are barred / total number of galaxies in that type). Errors are 
95% Bayesian binomial confidence intervals ||Cameron|2011[|. 


(Csabai et ai. 20031. Steiiar veiocity dispersions are taken 


(20111. An excerpt of these data may be 


from Oh et al. 
found in Table 


3.1 Barred AGN fraction at a fixed mass and 
colour 

Figurej^shows the distributions of mass and colour for AGN 
and star-forming activity types, split into barred and un- 


ence is primarily caused by selection effects relating to the 
underlying Eddington ratio distribution. The probability of 
a galaxy hosting an AGN is assumed to be independent of 
stellar mass, and thus AGN are prevalent at all masses in 
the range 9.5 < log(M/M0) < 12, despite only being observ¬ 
able at higher masses. As a result, we expect higher absolute 
numbers of barred AGN in a flux-limited sample since barred 
disc galaxies are also on average redder and more massive 
than unbarred disc galaxies ( Masters et al.|2011 |2012| ). We 
interpret this as the primary cause for the higher fraction 
of barred AGN (51.8%) versus barred star-forming (37.1%) 
galaxies in Table 

To control for this selection effect, we examine the 
fraction of AGN at fixed masses and colours (Figure [^. 
The total disc galaxy sample spans a mass range from 
9.0 < log(M/M0) < 11.5, while the colour range extends 
from 1.0 < (u — r) < 3.5. AGN hosts are found throughout 
the disc galaxy sample, but most appear in galaxies with 
log(M) > 10^° A70. When examining the fraction of galax¬ 
ies with an AGN as a function of mass and colour, redder 
and more massive galaxies have AGN fractions as high as 
10%. Bins with fewer than 10 total AGN (barred AGN -|- 


© 0000 RAS, MNRAS 000, 000-000 


















































Galaxy Zoo: bars and AGN 7 



Figure 4. Mass and colour distributions for disc galaxies in the GZ2 sample, separated by both activity type (either AGN or star¬ 
forming as in Table 1^ and the presence of a galactic bar. AGN (green) are on average both significantly redder and more massive than 
star-forming galaxies (blue). When splitting the disc galaxies into barred (solid lines) and unbarred (dashed lines), however, there is no 
significant difference between the two populations. Gounts are normalized so that the sum of bins is equal to 1 for each sample. 


unbarred AGN) are masked to minimize variance from small 
sample sizes. The same trend is also seen when splitting the 
disc galaxy sample into barred and unbarred subsamples. 

To analyze the difference between the barred and un¬ 
barred AGN popnlations, we plot the difference in barred 
and unbarred AGN fractions in Figure This quantity is 
defined as: 


du-NB = barred AGN fraction — unbarred AGN fraction 

( 1 ) 

and is calculated in each of the mass/colour bins in Figure]^ 
For each bin, a positive value represents a greater fraction of 
barred AGN and is coloured blue; a negative value represents 
a greater fraction of unbarred AGN and is coloured red. 

Since our AGN sample is divided into relatively small 
subsamples, we examine how the size and placement of the 
mass/colour bins affect the results of Figure]^ To control 
for this effect, we examine the average value of de-NB and 
the fraction of bins with dn-NB > 0, defined as: 


number of bins with higher barred AGN fraction 
total number of bins 

( 2 ) 

We compute /b>nb for 400 combinations of mass and colour 
bin widths between 0.2 Alog(M/M0) ^ 0.6 and 0.12 ^ 
A{u — r) ^ 0.35. The distribution of results from all combi¬ 
nations is shown in Figure]^ Our final bin choice (as seen 
in Figure]^ has a mass width of Alog(M/M0) = 0.375 
(16 bins) and colour width of A(u — r) — 0.16 (22 bins). 
This choice lies near the peak of the distributions for both 


/b>nb and du-NB, while maximizing the total number of 
bins to decrease the uncertainty on statistical tests. 

For the first time among recently published studies, we 
quantify the level of correlation between the presence of a 
bar and AGN through statistical analysis. We test the null 
hypothesis that in the absence of a causal link, the difference 
between barred and unbarred AGN fractions when binned 
by mass and colour should be centered around zero. The 
null hypothesis also requires that the likelihood distribution 
decreases symmetrically from zero in both directions; as a 
result, we assume a normal distribution with mean fi — 0 
and standard deviation a. Other models of the null hypoth¬ 
esis could of course also be tested, but we adopt this as the 
simplest reasonable scenario that fits the constraints of the 
problem. 


To assess the level of statistical significance, we fit 
the data in Figure with a range of models with vary¬ 
ing mean (^b-nb) and standard deviation (ad) and then 
apply an Anderson-Darling test. We selected this test be¬ 
cause it has been empirically shown to be more power¬ 
ful and reliable at testing normality than traditional or 
Kolmogorov-Smirnov tests, especially with small (n < 30) 
sample sizes jHou et al. 20091. The confidence threshold 
required for the model to pass at fitting the data is 95%. 
In Figure we show the distribution of the Anderson- 
Darling statistic as a function of ad for two of the tested 
models: the null hypothesis (du-NB = 0) and the best fit 
to the data (dn-NB = 0.012). The null hypothesis fails 
the Anderson-Darling test for all values of ad, indicating 
that the 66.7% ^216% fraction of bins that have a higher 
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Figure 5. Optical colour vs. stellar mass for disc galaxies in GZ2. Black contours represent all disc galaxies (top), all barred galaxies 
(middle), or all unbarred galaxies (bottom). All AGN (top), barred AGN (middle), and unbarred AGN (bottom) are plotted in the left 
panels as blue dots; the right panels show the AGN fraction in each colour/mass bin. Bins with A^/iGiV < 10 are masked. 


barred than unbarred AGN fraction is statistically signifi¬ 
cant. The best fit to the data, by contrast, has a mean of 
de-NB = 0.012^°'°°!^ and ad = 0.028. The positive value 
of de-NB indicates an increase in the AGN fraction for 
barred galaxies, consistent with the hypothesis that at least 
some fraction of AGN activity is triggered or sustained by 
bar-driven fueling. 


3.2 Comparing barred and unbarred AGN 
accretion strengths 

If the presence of a bar does contribute to AGN fueling, 
one possible result would be an increase in the accretion 
rate for barred AGN hosts vs. those that are unbarred. To 
assess this, we compare relative accretion strengths using 
the quantity R — L[o iiq/AfsH, with L[o iii] as a proxy for 
the AGN bolometric luminosity. [O ill] luminosities were 


calculated using fluxes from Oh et al. (20111, and black hole 
masses estimated using the Mbh-u relation: 


log 


Mbu 5 


,31og( 


200 km s' 




(3) 


Here a and j3 are empirical values determined from the ob¬ 
served relationship between black hole mass and velocity dis¬ 
persion a. We adopt the parameters measured by |Gultekin| 


et al. (20091 of {a,l3) = (8.12 ± 0.08,4.24 ± 0.41). 


It has been demonstrated for smaller samples of galax¬ 
ies that the parameters a and jS vary as a function of mor¬ 


phological type (Graham et al. 2011 Giiltekin et al. 2009 
[Brown et al.|2013| ), including differences between barred and 
unbarred galaxies. We choose not to use (a, /3) parameters 
where (a, P) are derived from separate subsamples for two 
reasons. First, since the Msu-a relation is calibrated from 
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Figure 6. Optical colour vs. stellar mass for barred and unbarred disc galaxies in GZ2. Coloured bins show the difference between 
the AGN fractions for barred and unbarred galaxies. Blue bins have higher fractions of barred galaxies, red bins have more unbarred 
galaxies, and pale/white indicates no difference. The region on the colourbar enclosed by the dotted lines represents the mean of the 
data determined by the Anderson-Darling test. The colour gradient is on the same scale as Figure]^ Bins with Ny^Qpf < 10 are masked. 
A colour version of this plot may be found in the electronic edition of the journal. 


small samples of nearby galaxies, the statistical error on the 
parameters increases as galaxies are divided into smaller 


sub-groups. The calibration of Giiltekin et al. (20091, for 


instance, is based on measurements of only eight barred 
galaxies. The error in /3 for the barred Mbh-ct relation is 
cr/3 = ±0.751, almost twice the error obtained by fitting to 
the full sample of disc galaxies. Second, while different stud¬ 
ies report consistent values for a and /3 when all disc galaxies 
are considered, the values can vary significantly when split¬ 


ting by morphological type. Lee et al.| (2012 ( and Alonso 
|et al.| ( [2013[ ) use separate values for (q,/ 3) and report con¬ 
flicting levels of agreement, depending on which parameters 
are used. This raises the possibility that differences in AGN 
strength are simply due to differences in calibration parame¬ 
ters, and not in the true distribution of accretion efficiencies. 

Figure|^shows the relative accretion strengths R for our 
sample as a function of mass and colour for both barred and 
unbarred AGN; these values are inversely correlated with 
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Figure 8. Fits of the binned fraction of barred vs. unbarred AGN fractions to a normal distribution. Left: value of the Anderson-Darling 
test (A^) as a function of the standard deviation of the normal distribution being fit (cti). The horizontal black line shows the critical value 
of corresponding to 95%; a model must fall below this line to be considered an acceptable fit at this level of confidence. Two models 
are shown: the null hypothesis (blue diamonds) and the best fit to the data in Figure]^ (purple triangles). Right: Plot of the minimum 
for the full range of means (cIb—nb) tested for the data. This shows that acceptable fits can be found for 0.005 <dB-NB< 0.019, but 
that the null hypothesis is ruled out at 95% confidence. 


both mass and {u — r) colour. This trend is likely driven 


by the same selection effects described in iS.l (Aird et al. 


20121. At a fixed L[o iii]/Mbh ratio, AGN with lower mass 


black holes are less likely to be detected due to the signal 
to noise requirements on their spectral lines. This biases the 
distribution of R toward higher mass black holes. Since stel¬ 


lar mass is strongly correlated with black hole mass (Haring 


fc RixpOOd Giiltekin et al.|[2009 Merloni et al.||2010 |, and 


stellar mass correlates with optical colour (Kauffmann et al. 


2003bI, this explaius the trend seen in both parameters for 


an uncorrected sample. 

Since these observationally-driven selection effects are 
likely to affect barred and unbarred galaxies equally, we 
compare the values of R of both groups without any cor¬ 
rections. A two-sided KS-test yields a p-value of p = 0.127 
for the two distributions. This is consistent with both the 
barred and unbarred galaxies being drawn from the same 
distribution. We thus conclude that there is no strong evi¬ 
dence for a difference in accretion strength between barred 
and unbarred AGN. 


This result contradicts Alonso et al. (20131, who found 


an excess of barred AGN with high values of R. We conjec¬ 
ture that this may be the result of their sample selection, 
which excluded galaxies with M* < lO^^M© in favor of a 
higher redshift limit of a = 0.1. However, low mass galax¬ 
ies have higher Lyo hij/Mbh ratios and are more likely to 
be unbarred than their higher mass counterparts (Lee et al. 
|2012[ ). If this effect is real, it appears to be limited to high- 
mass galaxies (which themselves are subject to selection ef¬ 
fects due to the methods used to measure R). Additionally, 
[Alonso et al.|(|^13[) include composites and LINERs in their 


sample of AGN. If the activity from these galaxies is not pri¬ 
marily from black hole accretion, R is not a true proxy for 
accretion strength, and comparisons between barred and un¬ 
barred galaxies do not accurately probe differences between 
the two populations. To test this, we compare R distribu¬ 
tions for barred and unbarred composite -|- AGN -|- LINER 
galaxies with M* > 10^^Mq. Eor these galaxies, the differ¬ 
ence in the average values of R for the barred and unbarred 
samples is 0.09 {Lq/Mq)~^ (compared to a difference of 
0.06 {Lq/Mq)~^ when considering only AGN with no cut 
on stellar mass), and a KS-test for the distributions yields 
a p-value < 0.01, which agrees with the results of j Alons^ 


et al. (20131. We note that our results are consistent with 


Lee et al. (20121, who have a similar mass range to our sam¬ 


ple of disc galaxies, and do not include composites in their 
sample. 


4 DISCUSSION 

We have compared a sample of 353 barred Seyfert AGNs 
to 328 unbarred Seyferts and measure the potential corre¬ 
lation between the presence of the bar and the AGN. We 
find that at fixed mass and colour, AGN hosts show a small 
increase in the fraction of galaxies that are barred. The av¬ 
erage difference is dn-NB = 0.012, or roughly 16.0% of the 
average barred AGN fraction. We find no difference in the 
^[o iiij/AIbh ratio between barred and unbarred AGN at 
either fixed mass or colour. We conclude that while AGN 
hosts have moderately higher probabilities of hosting a bar. 
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Figure 9. Left: Relative accretion strength R vs stellar mass for barred (blue) and unbarred (red) AGN in our sample. R is plotted as 
the mean of values within five equal-width bins in the range 9.8 < log(M/M 0 ) < 11.3, which includes 98% of the AGN sample. Points 
are drawn at the midpoint of each bin. Right: R vs colour for barred and unbarred AGN. R is plotted as the mean of values within 
five equal width bins in the colour range 1.6 < {u — r) < 3.0, which includes 96% of the AGN sample. Error bars for each plot are 95% 
confidence intervals, calculated by bootstrapping with 1000 times resampling. There is no significant difference in accretion strengths for 
barred and unbarred AGN as a function of either mass or colour. 


the presence of the bar does not seem to affect either the 
quantity or efficiency of fueling the central black hole. 


If bars are not required to initiate AGN fueling, then 
what is the source? There must be a process that trans¬ 
ports angular momentum through the galactic disc and 
creates/maintains an accretion disc. Both theoretical mod¬ 
els i Shlosman et al.||1989 19901 and numerical simulations 
(Hopkins & Quataert 20101 indicate that this process re¬ 
quires two stages. First, the gas must be driven from a 
radial scale of megaparsecs down to kiloparsecs. Standard 
viscous torques on the gas are too inefficient to initiate gas 


inflow by themselves Shlosman et al. (19891; Bournaud et al. 


120051; therefore, some other mechanism is required. Within 
the central kiloparsec, a secondary instability must take over 
within the gaseous disc for AGN fueling to occur. 


In the context of this general model, we consider three 
possibilities: (I) Bars are a necessary ingredient for fueling 
AGN, (II) Bars are one of multiple processes that fuel AGN, 
or (III) Bars play no role in fueling AGN. We also discuss in 
each scenario possible explanations for the existence of all 


four observed combinations: barred AGN, unbarred AGN, 
barred non-AGN, and unbarred non-AGN. 


4.1 Scenario I: Bars are necessary to fuel AGN 


If the presence of a stellar bar is the only mechanism by 
which gas can be driven to the ~ 1 kpc scale, there must 
be a reason both barred and unbarred AGN are observed in 
large numbers. One possibility is that a galactic bar initiates 
fueling of the black hole, but is subsequently destroyed in 
a dynamic timescale shorter than the lifetime of the AGN. 
These separate timescales are not currently known with cer¬ 
tainty, but estim ates place the lifetime of an AGN from 10® 
— 10® years (eg, Schawinski et al.|2010 Martini|2004 1. The 
range of bar lifetimes is not yet firmly established; some 
models show bars to be transient features that are destroyed 
either due to buckling from angular momentum transport 
or from the build-up of a central mass concentration (CMG) 
(Bournaud et al.|2005 Combes||2008 1. In these models, the 
lifetime of a bar is estimated to be 1 — 2 Gyr. [Kraljic et al.| 
(|2012[) also found bars to be short-lived in their simulations. 
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Figure 7. Distributions of the difference in the fraction of bins 
with excesses of barred AGN (/b>nb) the average dif¬ 

ference between barred and unbarred AGN fractions (dB-NB)- 
Both values are computed for 400 variations in the mass and 
colour bin widths. Left-. The average fraction of bins with a 
higher barred AGN fraction is /b>nb = 0.705 it 0.073. Right: 
The average difference in barred and unbarred AGN fractions is 
<^b—NB = 0.015 it 0.004. Dashed black lines indicate the values 
of /b>nb and average dB—NB used in Figure and subsequent 
analysis. 


4.2 Scenario II: Bars are one of several ways to 
fuel AGN 


If stellar bars are only one of several ways to fuel AGN, 
then both barred and unbarred AGN should exist (as should 
both barred and unbarred star-forming galaxies). The sim¬ 
ulations conducted by Hopkins & Quataert ( 2010| support 
this model, which show that multiple large-scale mechanisms 
(including a stellar bar) can be responsible for transporting 
gas to scales required for AGN fueling. Further, if bar-driven 
fueling is responsible for some fraction of the AGN, this 
model predicts an increase in the fraction of barred AGN, 
which our data supports. 

While the existence of unbarred AGN is explained by 
this model, there is no immediate explanation for the ex¬ 
istence of barred galaxies that do not host AGN; here we 
suggest several possibilities. First, a bar that initiates AGN 


fueling may simply outlive the AGN (see 4.11, which agrees 
with estimates of both bar and AGN lifetimes. Second, there 
could be a correlation between bar strength and AGN activ¬ 
ity, where only sufficiently strong bars initiate fueling. This 
is consistent with Lee et al. (20121, who find a higher AGN 


fraction in barred galaxies where the bar length is at least 
1/4 of the total disc diameter. They did not test, however, 
whether this relationship remains at fixed mass and colour. 
Finally, the emission from an AGN is expected to be highly 
variable with time, driven by processes such as accretion disc 
instabilities and/or feedback within the accreting material 
( Hickox et al.|[2bl4| . In this case, barred galaxies without 
AGN are simply observed in low parts of their duty cycle, 
with Eddington ratios too low to be detected at the limits 
of our observations. 


4.3 Scenario III: Bars do not fuel AGN 


but only early bars (formed at z > 1). Bars formed later (at 
z < 1) were maintained down to z = Q, giving a lifetime of 


Finally, we consider the possibility that stellar bars do not 
trigger AGN activity in any way. This is inconsistent (al¬ 
though marginally so) with the increase in barred vs. un¬ 
barred AGN fractions that we find at fixed mass and colour. 
One possibility is that the model used for the null hypoth- 


at least 8 Gyr. 



esis (a normal distribution centered at dB-NB= 0) does not 

Other simulations ( 

Debattista 

et al. 2004 2006 

apply. Detailed simulations of cosmological volumes that in- 

Athanassoula et al.||2005 

2013 Shen 

& Sellwood||2004 1 do 

elude both AGN and detailed disc morphology, such as II- 


not observe bar destruction due to buckling. In these cases, 
only a sufficently massive CMC is capable of destroying the 
bar on Gyr timescales. The mass of the CMC required in 
these models is at least several percent of the total mass of 
the disc — this is significantly larger than the mass measured 
in local disc galaxies. If the CMC is insufficently large, the 
bar is maintained for the lifetime of the disc (up to 10 Gyr; 
Athanassoula et al.|2013 1, and thus should be observable for 


at least the lifetime of the AGN. 

If bars are truly long-lived structures in all disc galaxies 
and are necessary to fuel AGN, we would expect a much 
higher value of the ratio of barred to unbarred AGN hosts. 
Since the observed numbers are nearly 1:1, we consider this 
scenario highly unlikely. It is possible that bars are necessary 
to fuel AGN, but the number of observed unbarred AGN can 
only be explained if the factor of ~ 10 difference between the 
upper end of the AGN lifetime and the lower end of the bar 
lifetime can be resolved. While this is possible, we consider 
it unlikely given the assumptions required. 


20141 should ultimately provide more well-defined priors for 


this. 

In addition, our test of the null hypothesis could still be 
consistent with a strong effect even if the total number of 
barred and unbarred bins were equal. For example, if bar- 
driven fueling is strongly mass-dependent, the ^b-nb bins 
could have excesses of barred AGN at high masses and 
deficits at low masses; this would still be consistent with 
a distribution centered at zero. We test the simplest cases 
by simply splitting the sample into two in both mass and 
colour (Tablej^. Low- and high-mass disc galaxies (dividing 
the sample at \og{M/M q) = 10.625) have nearly identical 
values of /b>nb and mean ^b-nb; there is no evidence of 
a mass-dependent effect on bar-driven AGN fueling. When 
splitting discs into red vs. blue (at a colour of (u—r) = 2.22), 
bluer galaxies do have significantly more bins with an ex¬ 
cess of barred AGN (/b>nb= 0.88) than redder galaxies 
(/b>nb= 0.54). The uncertainties on /b>nb are quite large. 
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sample 

/b>nb 

Mean c^b—NB 

low mass 

log(M/M 0 ) < 10.625 

0.70 

0.0125 

high mass 

log(M/M©) > 10.625 

0.64 

0.0123 

blue 

(u — r) < 2.22 

0.88 

0.023 

red 

(u — r) > 2.22 

0.54 

0.006 


Table 3. Difference between barred and unbarred AGN fractions 
for disc galaxies when splitting the sample in two by both mass 
and colour. /b>nb is the fraction of bins that show an excess of 
barred AGN (compared to unbarred), while de-NB is th® average 
value of the differences over all bins. Since the number of bins in 
each subsample is only ~ 8 — 13 when splitting by mass or colour, 
the uncertainty in /b>nb is correspondingly large. 


though, since each subsample has less than a dozen bins. 
Our splits by colour agree with ( Oh et aL]|2012 1, who find 
that bar effects on AGN are more pronounced in bluer and 
less massive galaxies. |Lee et cJ!] ( | 2012[ ), in contrast, find that 
/b>nb depends on neither mass nor colour. 

If bars have no impact at all on the likelihood of a disc 
galaxy hosting an observable AGN, this is inconsistent with 
both the models and simulations that demonstrate efficient 


gas-driven inflow by bar structures (Hopkins & Quataert 
|2010[ ). If the efficiencies of other morphologies that drive 
gas inflow are much higher than bars, though, this could 
also be consistent with our data. A lack of bar-driven fuel¬ 
ing is consistent with the existence of both barred and un¬ 
barred AGN and star-forming galaxies, and the nearly equal 
numbers found in both pairs. 

Given the limits on the data set (which is driven by bin¬ 
ning the total number of disc galaxies by mass and colour), 
we do not completely rule out this model. However, given 
the small (but measurable) increase in the bar fraction from 
our data and the current constraints on both bar and AGN 
timescales, we propose that bar-driven fueling must account 
for at least some fraction of observed AGN activity (^4.2[|. 


5 CONCLUSIONS 

We have created a sample of 19,756 disc galaxies from SDSS 
DR7, using data from the Galaxy Zoo 2 project for morpho¬ 
logical classifications of strong, large-scale bars. We studied 
the effects of stellar bars on 681 AGN and compared these 
effects to a control sample of disc galaxies both without bars 
and without AGN. The Galaxy Zoo 2 data provides a very 
large sample of disc morphologies for which the bar likeli¬ 
hood can be empirically quantified, based on crowdsourced 
visual classifications. 

We find that the fraction of barred AGN (51%) is sig¬ 
nificantly greater than the fraction of barred galaxies with 
central star formation (37%). However, this is driven both by 
selection effects for detecting optically-identified AGN and 
by known correlations between black-hole mass and stellar 
mass, as well as stellar mass and optical colour. When ex¬ 
amining the fraction of barred AGN as a function of a fixed 
mass and colour, we still find a small increase in the number 
of barred AGN hosts. The null hypothesis of no relationship 
between the two cannot be ruled out at the 95% confidence 
level. The L[o iii]/Mbh ratio R (a proxy for the overall ac¬ 


cretion rate) shows no dependence on the presence of a bar, 
once the same mass and colour constraints are applied. 

Our results are consistent with a small relationship be¬ 
tween the presence of a large-scale galactic bar and the pres¬ 
ence of an AGN. We propose that while bar-driven fueling 
does indeed contribute to some fraction of the current ob¬ 
served population of growing black holes, other dynamical 
mechanisms, such as lopsided or eccentric stellar disk, must 
also contribute to the redistribution of angular momentum 
and thus the fueling of the accretion disk at small galactic 
radii. 

Even with the advent of the large-scale SDSS data and 
the morphological classifications from Galaxy Zoo 2, this re¬ 
sult is still constrained by the total number of galaxies in our 
study. Larger samples of disk galaxies with activity and mor¬ 
phological classifications, notably the Dark Energy Survey 
(DES) and the Large Synoptic Survey Telescope (LSST), 
should increase the sample sizes by factors of at least a few 
and help to confirm these results. Further development on 
the theoretical side is also critical — with state-of-the-art 
simulations now able to reproduce both the morphology dis¬ 
tributions and the observed black hole mass function, these 
results can be compared to theory in a cosmological context. 
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